FOURIER DUALITY FOR FRACTAL MEASURES WITH AFFINE SCALES 



DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 



Abstract. For a family of fractal measures, we find an explicit Fourier duality. The measures in 
the pair have compact support in R d , and they both have the same matrix scaling. But the two use 
different translation vectors, one by a subset B in R d , and the other by a related subset L. Among 
other things, we show that there is then a pair of infinite discrete sets T(L) and T(B) in R d such 
that the r(L)-Fourier exponentials are orthogonal in L 2 (hb), and the F(_B)-Fourier exponentials 
are orthogonal in L 2 (/il). These sets of orthogonal "frequencies" are typically lacunary, and they 
will be obtained by scaling in the large. The nature of our duality is explored below both in higher 
dimensions and for examples on the real line. 

Our duality pairs do not always yield orthonormal Fourier bases in the respective L 2 (/i)-Hilbert 
spaces, but depending on the geometry of certain finite orbits, we show that they do in some cases. 
We further show that there are new and surprising scaling symmetries of relevance for the ergodic 
theory of these affine fractal measures. 
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1. Introduction 

Fractal scaling and self-similarity occurs in nature, and in applications, such as to large com- 
munication networks. To understand the nature of fractal scaling, it has proved useful to develop 
specific model cases. Here we focus on such a class, those given by a finite family of affine trans- 
formations. Iteration in the small yields fractal measures [/,, and their support sets fall into one of 
the classes of compact fractals embedded in M. d . Iteration of scale in the large, in turn, leads to 
a kind of fractal networks comprising of a lattice skeleton and lacunarity, degrees of sparsity. We 
study when these scales in the large lead to orthogonal families in L 2 ([i). 



Supported in part by the National Science Foundation. 

2000 Mathematics Subject Classification. 47B32, 42B05, 28A35, 26A33, 62L20. 

Key words and phrases. Affine fractal, Cantor set, Cantor measure, iterated function system, Hilbert space, 
Beurling density, Fourier bases. 

1 



2 



DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 



There have been a number of recent papers dealing with and treating a variety of features of the 
harmonic analysis of fractal measures \x with affine scales, see for example [DJ06a, DJ07a] IDJ07bt 
ID J07cl IDJOM IDJ081 IDJ091 IGIL091 IHL081 HW081 IStrM ISW991 IStrOO] . 

In this paper we study pairs of orthogonal complex exponentials in two Hilbert spaces L 2 (n) for 
a pair of fractal measures [i with affine scales. Each \x in the pair has the same scaling matrix, but 
the affine mappings making up the iterated function systems (IFSs) are different. 

The particular pairs we study are selected from a certain axiom involving a fixed complex 
Hadamard matrix. It is indexed by a triple, a matrix and two sets of vectors: more specifically, by 
an expansive d by d matrix R with integer entries, and by two subsets B and L in M. d having the 
same cardinality. We show that when the data (R, B, L) are fixed, we then get a naturally defined 
pair of affine measures \ib and hl, each with i?-selfsimilarity. The two measures arise by taking 
scaling in the small with powers of the inverse R , and initiating with the given sets B and L. The 
support of the measures is typically a Cantor fractal, e.g., a Cantor set on the line, or for example 
a planar Sierpinski set in M 2 . 

We offer a detailed Fourier duality for the pair, and we show that there is a pair of discrete sets 
r(L) and T(B) in M. d such that the r(L)-Fourier exponentials are orthogonal in L 2 (//#), and the 
r(£?)-Fourier exponentials are orthogonal in L 2 (iil). These sets of orthogonal "frequencies" will 
be obtained by scaling in the large. The nature of this duality is explored below both in higher 
dimensions and for examples on the real line. 

Our duality pairs do not always yield orthonormal Fourier bases in the respective L 2 (/i)-Hilbert 
spaces, but we show that they do in some cases. We further show that there are new and surprising 
scaling symmetries of relevance for the ergodic theory of these affine fractal measures. 

The paper is organized as follows. The problem we consider here began with the question 
raised first by Fuglede |Fug74 for open subsets O in M. d concerning orthogonal Fourier bases for 
the corresponding Hilbert space L 2 (Q) with respect to Lebesgue measure, i.e., of orthogonal bases 
(ONBs) of complex exponentials. While the possibility of such Fourier bases is related to tiling 
properties for f2, it is known not to be equivalent [Tao04]. For a brief summary of these questions, 
see for example [JP991 |IK P99J |ffi T03 j EZPgH [J0r82b| [J0r82a|. 

However, the ideas from [JP99J and [DJ07a] still suggested intriguing connections between spectra 
and geometry, more specifically that Fourier bases should be tied to an intrinsic selfsimilarity, and 
further that this can be made precise with the use of Hadamard matrices (see Definition [23]). But in 
addition, this suggested algorithmic iterations of the similarity transformations, so the introduction 
of iterated function systems (IFSs) built from a finite family of affine transformations in M. d (see 
Definition I2.ip . Thus, the measures /i arising in the limit can be expected to come from a dual 
affine iteration. Further [JP99] suggested a definite relationship between the two sides in such a 
duality. 

In section [2] we introduce the Hadamard matrices, and we prove the stated relations for the two 
IFS-measures in duality. In Proposition 12.61 we show that each of the Hilbert spaces L 2 (fi) for the 
dual system of measures have a natural pair of infinite families of orthogonal Fourier frequencies. 
The fractal measures arise by iteration in the small, while the orthogonal families by iteration in 
the large. 

The interplay between the two sides in the duality is made precise in Corollaries 12.81 and 12.91 In 
Theorem 12.101 we introduce a class of invertible matrices which define equivalence between the two 
measures. Further we offer examples in section [3] to the effect that the two Hilbert spaces L 2 (/j,) 
for dual measures may have different Fourier properties. 
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Further we show in Lemma [3 .31 that the nature of the orthogonal families from Proposition 12.61 is 
determined by geometry. By this we mean that for infinite families of orthogonal Fourier frequencies, 
the obstruction to the ONB property is determined by whether a certain pair of dynamical systems 
have non-trivial orbits. A further aspect of the duality is pointed out in Lemma 14.11 with the use 
of a pair of transfer operators. 

We offer computable examples in section [3] and [5] which we believe may be of independent 
interest. The example in section [5] has already been the subject of extensive work, see for example 
[JP98llDJ09llDHS09j . 

In section [5] we prove that L 2 (n) may have orthogonal Fourier bases ONBs (i.e., ONBs of complex 
exponential frequencies) having arbitrarily small upper Beurling density. Compared with what is 
usually expected in standard sample theory, this result reveals new features of the harmonic analysis 
of the Hilbert spaces L 2 ([i). 

2. Pairs of spectral fractal measures 

There is a general construction of measures with intrinsic selfsimilarity due to Hutchinson 
[HutST], but it is closely related to families of infinite product measures considered earlier. The 
starting point for this construction is a finite family F of contractive transformations in a complete 
metric space and an assignment of probabilities on F. Repeated iteration of the mappings in F, and 
taking averages, then leads to a unique Borel probability measure \i in the limit; see f)2.2f) below. A 
special case of this is the measure resulting from a Cantor iteration, and recursive rescaling. Mo- 
tivated by this, we will be concerned here with the special case when the mappings in F are affine 
transformations in M. d for some d, see (j2.4|) . In that case, there is a family of complex exponentials 
indexed by points in M d , and we will be interested in iterative algorithms for the construction of 
maximal orthogonal families in L 2 (/i). Ideally we ask for these families to form orthonormal bases 
(ONBs) for L 2 (/i), i.e., Fourier bases. 

Definition 2.1. Let (Y,d) be a complete metric space, and let (ri)fL i be a finite system of con- 
tractive mappings, i.e., there is a constant c, < c < 1 such that 

(2.1) d(Ti(x),Ti(y)) <cd(x,y), (x,y € Y,i = 1, . . . ,N) 

Let {pi)f-i, Pi > satisfy Y2i=i Pi = 1- Then by |Hut81] there is a unique probability Borel measure 
\i such that 

N 

(2.2) f = ^w oT r 1 

i=l 

where \x o r^ 1 is the measure given by (/x o t~ 1 )(E) = fj,(T~ 1 (E)), E G Q3(y)=Borel sets, and 
r -i(£) = { x G Y : t(x) eE}. 

Here we will study the case Y = W 1 with its usual metric, and we will be interested in families 
(tj) consisting of transformations of the form 

(2.3) x ^ R~\x + b) 

where R is a fixed expansive matrix (i.e., all eigenvalues A satisfy |A| > 1) and where b is in a finite 
subset B of R d . 

We will further restrict to the case of equal probabilities pi = 1/N in (|2,2p . 
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The mappings in (|2.3p will be denoted 

(2.4) r b (x) = R- 1 (x + b) 

It is known that the corresponding measure \i = hr,b solving (|2.2p has as support a Cantor set 
R d . However, depending on B, the measure \x will vary. It is convenient to chose B such that 
£ 6. In that case 

(2.5) supp(//) = \jr R-% : b k G b] =: X(B) 
Setting 



e t (x) := e 2nit - x , (t,x£R d ) 

we get 

Definition 2.2. Let r C R d be some discrete subset, and let 

E(T) := {e 7 : 7 G T} . 
We say that F is orthogonal in L 2 (fi) iff the functions in E(T) are orthogonal, i.e. 

(e 71 , e 72 ) = £(72 - 71) = 0, for all 7l / 72 G T, 
where ju is defined as the Fourier transform 

(2.6) At(*) = / et(x)dn(x). 



Definition 2.3. Set N := #B and 

(2-7) XB (t) = l^ efe (t) 



N 

beB 

(2.8) (5b := — ^2$b (Dirac notation); 

beB 

then 

(2.9) XB(t) = 5 B {t) 

Lemma 2.4. Let R and B be as above, and let (t^^b be the IFS in (|2.4p . Let jj, = Urb be the 

corresponding Hutchinson measure, and let R T be the transposed matrix. Then 

oo 

(2.10) m = U^{(R T r k (t)] 

k=l 

where the infinite product is absolutely convergent. 

Proof. Well known! □ 

Definition 2.5. Let a d x d matrix R be given. Assume R G M^(Z) and that R is expansive. Let 
S, L C M d be such that G S, G L, JV = #5 = #L and assume 

(2.11) R k b-le Z, for all 6 G B, / G L, k G Z, fc > 
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and further that the matrix 

1 f2mR- 1 b-l 



(2.12) (e 2 , 

is unitary. 
Set 

(2.13) T(B) ■= \J2 Rkbk -h^ B,neZ + \ , and 

U=o J 

(2.14) r(L) := I ^(^ T )^fc : ifc e L, n G Z + 1 

U=o J 

Proposition 2.6. Let be the Hutchinson measure for (Tf,)beB and let I^L be the Hutchinson 
measure for the dual system 

(2.15) o-i(x) := (R T )-\x + l), (I G L) 
Then T(B) is orthogonal in L 2 (fii), and T(L) is orthogonal in L 2 (/x#). 
Proof. The conclusions of the proposition amount to the following: 

(2.16) (e 7l , e j2 )^ B = fi B (li - 72) = £71,72 for a11 7l,72 G T(L) 
and 

(2-17) (e 6 , e 6 ) Mi = - 6) = f ° r a11 &>fc G r(S) 

In view of (12. lip and (|2.12p in the definition, it is easy to prove one of the stated properties and 
we offer the details in the verification of (|2.17p . Using Lemma 12.44 we get 

00 

(2-18) MZ) = HxL(R- k Z) 

k=l 

and we will be using this for £ = £1 — £2 / 0, with pairs of points from the set T(B). So we are 
concerned with 



^ = Y^R\b k -h 



fc=0 

where b , . . . , b n , f3 , . (3 n G B. 

Substitution into (|2.18p shows that the product representation for — £2) will contain a 

factor of the form 

(2.19) xdR-\b-/3)) 
where b,(3 G B distinct. Then 

(2.20) XL(R-\b -/?)) = 1 J2e^- 1 ^e 2mR - lM = (row^ , row 6 ) 

leL 

where row^ is a notation for the row with index j3 in the Hadamard matrix (|2.12p . and where (• , •) 
in (|2.17p refers to the standard complex inner product in C . 

Since the matrix in (|2,12p is assumed unitary, it follows that distinct rows are orthogonal. The 
desired conclusion follows. □ 
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Corollary 2.7. Let the data (R,B,L) be as stated in Proposition \2. (k then both of the orthogonal 
systems T(B) and T(L) are infinite. 



Proof. By symmetry, it is enough to prove this for T(B). We will show that if 60,61,... and 
Po, f3±, . . . in B yield the same vector 

n 

k=0 k=0 

then bk = (5k fo r < k < n. 

If not, then let k be the first term in (|2.2ip with bk 7^ Pk- Then, from (|2.2ip . we get 



(2.21) 1 = Y J R k b k = Y, Rk h 



R-\b k -h)= E R^-^i-h) 

i=k+l 

and, using (I2.1ip this implies that 

e 2,iR-H k 4 = e 2*iR-^ k ^ for aU l € L 

But then the rows bk and @k in the Hadamard matrix (|2.12p coincide, so they cannot be orthogonal. 
The contradiction implies the corollary. 

□ 



Corollary 2.8. Setting 

J 



(2-22) 4fl)(*)= E l^(t + 7)| 2 , 



(2-23) 4tt)(*) = E \Mt + 7)\ 2 , 



7er(L) 

and 

r {L) , 

?er(B) 

we note that the two functions are entire analytic oft = (t\, £2, . . . , to), i.e., have entire analytic 
extensions to C d . 
Further we have: 

o- { r f L) (t) < 1, and ffftL(t) < 1 for all t G R d . 

Proof. While most of the conclusions follow from earlier papers, we include here the proof for 

in (|2.22p . The analogous formula holds for o~^ B ^ in ()2.23p . 

Let Pl denote the orthogonal projection in L 2 (fj,s) onto the closed subspace 7i(L) spanned by 
E(T(L)), i.e., by the exponentials 

E(T(L)) :={e 7 : 7 er(L)}cL 2 (/iB). 

Then 

4?1)W= E \VB(t + 7 )\ 2 = E \^-t^~ f ), B \ 2 = \\PLe-t\\\ HllB) 
7G r(L) 7 er(L) 

= (e- t , P L e- t ) L 2 {llB) = (eo, U(t)P L U(-t)e ) L2(jlB) , 

where, for t = (t\, . . . ,td) £ U(t) denotes multiplication by e 2m ( tixi ~ { ^d^d). Since supp(^s) = 
X(B) is compact, U(£)PiXJ(—t) is entire analytic. 



□ 
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The computation further shows that 

4*i)(t) < \\U(t)PLU(-t)\\\\eo\\l Hw) = \\Pl\\ = 1 
where we used the fact that Pl is a projection so has norm 1. 

Corollary 2.9. The following conclusions hold: 

(i) T(L) is total in L 2 (fi B ) (i.e., it is an ONB) iff a^L) = 1 in Rd - 

(ii) T(B) is total in L 2 (fi L ) iff af) B) = 1 in R d . 

Proof. See [JP99] . [HTM] . □ 

Theorem 2.10. Let R, B, L be as stated in Definition \2.5[ Assume further that R T = R. If there 
is aGe GL d (R) such that G = G T , 

(2.24) GR = RG 

(2.25) G(B) = L 

then the two spectral functions f|2.22j) , (|2.23p in Corollary \2.8\ satisfy 

(2.26) GT(B) = T(L) 
and 

(2.27) (*) = a ?{L) ( Gt ) f° r aU 1 G Rd - 

Proof It follows from ||223jt and (j2T23jl that (12T26]) is satisfied. 
Since G(B) = L, we have 

x*<® = iyYl e o(6)(*) = ^ E e ^ Gt ) = ^ Gt ) 

fees befi 

and 

oo oo oo 

Mt) = n x L (R- k t) = n xB(Gii- fc t) = n x^-^) = ^go- 



k=l k=l k=l 

pd 



Furthermore, for t £ 

4(l)W= E \Mt + l)\ 2 = E lMB(G(t + 7 ))| 2 = E l^(Gt + G 7 )| 2 
7er(B) 7er(B) 7er(B) 

= £ |/2 B (Gt + 0| 2 = 4(L)(^) 
5er(i) 

which is the desired formula (j2.27|) . □ 
Corollary 2.11. If d = l, and #(B) = #(L) = 2, tten i/iere exists j£R \ {0} suc/t f/mi 

SL) ,.y _ JB) 



s 
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3. HADAMARD MATRICES AND EXTREME CYCLES 



In this section we prove that iteration by some fixed scaling matrix and its transposed leads to 
a new Fourier duality for a pair of IFS measures. Each of the two measures arises from an affine 
system, but the systems are linked by a complex Hadamard matrix, see Definition 12.51 and Example 
3. II below. Lemma |3 . 31 offers a geometric tool which allows us to test when our iterative algorithms 
from section [2] for the maximal orthogonal Fourier families in L 2 (fj.) in fact lead to orthonormal 
Fourier bases (ONBs) in the respective L 2 (^)-spaces. 

The basic building block for our examples is the matrix (|2.12p . An N x N matrix T~L over C is 
called Hadamard if \Hj,k\ = for all j, k = 1, . . . , JV. 

Complex Hadamard matrices serve as tools in combinatorics, in algebra, and in applications, 
see e.g. [CHK971 IDen09l IDit03l IDit04l IDFdGtHROU IDvA08l IPed041 |XCQ05| . The simplest are 
the unitary N x N matrices which define the Fourier transform on the finite cyclic groups Z^r = 
Z/NZ = {0, 1, . . . , N - 1}. If C is the iV-th root of 1, f = e 2ni / N , then 

(3-i) -)=(C J '-*W W 



is a complex Hadamard matrix. Moreover, the complex Hadamard matrices are closed under the 
following operations: 

(i) permutation of rows; 

(ii) permutation of columns; 

(iii) multiplication of a fixed row by a fixed phase; 

(iv) tensor product. 

By (iv) we mean this: if U is N x N and V is M x M complex Hadamard matrices, then 
W := U <g>V is also a NM x NM complex Hadamard matrix. To see this note that U V has 
entries of modulus l/y/NM. With the definition of^gFonC^® C M = 

(U ®V)(z<8> w) := Uz (g> Vw 



iNM 



we see that U C 
Example 3.1. 

then 



) V is again unitary. 
If 

V2 



U ' 



and V = 



1 

v 7 ^ 



By contrast, the Fourier transform of the group 



The main distinction between duality pairs \ib and /^l in the case d 
with the following: 
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Definition 3.2. Let T b {x) = R^ 1 (x + b) be an IFS as specified in (|2.4p . If X(B) is defined as in 
(|2.5p . a finite B-cycle point is a point in X(B) of the form (w,w, . . . ) obtained by the repetition 
of a fixed finite word w = (pi . . . b n ). When w is fixed, we denote by x(w) the corresponding finite 
cycle point, i.e., setting 

x n := R~ l b x + • • • + RT n b n 

we get 

x(w) =x n + R- n x n + R- 2n x n + .•• = (/- R-^Xn = R n (R n - Z)" 1 ^ 
Note that can also be defined by the property 

r bl ... n n x{w) = x(w) 

We call the set {x(w), Tf }n x(w), Tf )n _ 1 Tf )n x(w) , . . . , t& 2 . . . Tb n x(w)} the B-cycle generated by x(w). 

Similarly, one can define L-cycle points and L-cycles, associated to the IFS {(Ji)i & l (see (I2.15P ), 
Consider a system R,B,L in R rf satisfying the conditions in Definition 12.51 A finite B-cycle C in 
X(B) is said to be L-extreme if = 1 f° r an x £ C. A finite L-cycle C is said to be B-extreme 

if |xb(x)| = 1 for all x G C. 

Lemma 3.3. Lei R,B,L in R ('so £/ze dimension d = 1) be as specified in Definition \2.5[ In 
particular, we assume G i? and G L. 

(i) Then T(L) is an CWL in L 2 {^lb) iff the only B-extreme cycles in X(L) are the singleton 

!"}• 

(ii) Moreover, the set T(B) is an ONB in L 2 (fj,i) iff the only L-extreme cycles in X(B) are 
the singleton {0} 

For dimensions d > 1, the condition on the extreme cycles is only necessary, but not always suffi- 
cient for the corresponding T set to be an ONB. 

Proof |LW021 IDJ06b] □ 

In the analysis of extreme -B-cycles, the following lemma will be useful. See [DJ07a, Theorem 
4.1] for more details. 

Lemma 3.4. If C is a B-extreme cycle, then for all x G C, and all b G B, b ■ x G Z. 
Proof. For all x G C, we must have |x_b(^)| = 1- Therefore 



E 

beB 



N. 



All the terms in the sum have absolute value 1. There are of them, and one of them is 1, since 
G B. Therefore we have equality in the triangle inequality and this implies that e 2mb ' x = l for 
all beB. Therefore b ■ x G Z for all b G B. 

□ 



In the following, we point out the simplifications resulting from specialization to d = 1. In 
particular (Proposition 13. 5[) we point out that the dual systems from section [2] when specialized to 
the particular Hadamard matrices defining finite Fourier transforms yield intriguing pairs of fractal 
measures in duality. 
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In this section, we consider a family of examples which are related to the finite cyclic group Z^v, 
but they display fractal features and fractal duality which has not been studied earlier in duality 
theorems involving the finite cyclic groups. 

Setting 4.1. Fix integers M and N in Z+, and assume N\M, i.e., M = qN for some q G Z + , 
q > 1. We will consider the following instance of a system (R, B, L) in M. d subject to the conditions 
in Definition [231 Set R = M and L = {0, 1, . . . , N - 1}, and B = qL = {0, q, 2q, . . . , (N - l)q}. 

As a result, we see that the Hadamard matrix for R, B, L is 



(3.2) -4?fC W ) , C 



e 2iri/N 



Proposition 3.5. Let (R,B,L) be a system constructed from the numbers N and M as in Setting 
4-1. Then 

(i) T(L) is an ONB in L 2 {^ B ); 

(ii) T(B) is an ONB in L 2 (fi L ). 

Proof. We will establish (i) and (ii) as an application of Lemma 13.31 Note that the IFS which 
generates \xb is 

(3.3) A B \x) = X -±p-, keZ N ; 
and for fii it is 

(3.4) 7i W( x ) = £ + l > i e z N . 
For the two functions xb and XL we have 



iV-l 

(3-5) XL (t) = _^e fc (t) 

fc=0 

and 

(3-6) X B(t) = XL(qt). 

In exploring the extreme cycles we note that 

(3.7) |xL(t)| = lifft€Z; 
and 

(3.8) \x B (t)\ = imt€-Z. 

Q 

Since X(B) C [0, ^i], X(L) C [0, g^] and M = gJV, g > 1, it follows that the only L-cycle in 
X(B) is {0}. And similarly, the only i?-extreme cycle in X(L) is {0}. The result now follows from 
Lemma 13.31 □ 



The next example relates to the pair of fractals in Proposition 13.51 naturally associated with the 
finite Fourier transform of Z n . Since Z ra is its own Fourier dual, the corresponding B-L duality 
may be phrased in the language of the Z ra -Fourier transform, with the two sets B and L essentially 
being a copy of Z n . Example 13.61 below shows that if one of the sets is changed by one point, then 
the ONB conclusions from Proposition 13.51 no longer holds. 
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Example 3.6. The following example shows that if the numbers in the sets B or L in Proposition 
13.51 are modified, then the ONB conclusion may fail. 

Let R = 8 and B = 2 • {0,1,2,3} = {0,2,4,6} as in Proposition I3.5| but for L now choose 
L = {0, 1,2,7}, i.e., change 3 to 7 as compared to Proposition 13.51 with N = A and q = 2. Now 
x = 1 satisfies t^x = x and IxbO^)! = 1- Therefore {1} is a cycle in X(L) which is .B-extreme. 
Combining this with the earlier observations we see that the following hold: T(L) is not an ONB 
in L 2 ([i B ), but T(B) is an ONB in L 2 {[i L ). 

The next example shows how the two spectral functions in Corollary 12.81 will typically be quite 
different when the conditions in Theorem 12 . 101 are not satisfied. In the example the non-zero vectors 
in B are linearly independent while they are proportional in L. As a result there cannot be an 
invertible matrix G satisfying the conditions in Theorem 12.101 

Example 3.7. Let 




0.5 












/ 


0.4 














0.3 










/ 
/ 




0.2 






h. 

1. 1. 


/ 






0.1 






b L 
Ct. L.L 
















/ 

/ 
























0.1 0.2 


03 0.4 0.5 









FIGURE 1. The attractors X{B) and X(L) 

We prove that T(L) is a spectrum for but T{B) is not a spectrum for 
Set C := C3 = e 2ni/3 . Then the matrix in (I2TT2D is 

1 

C 

c 2 

which is the matrix of the Fourier transform on Z3. 

It is easy to see that X(B) C [0, |] x [0, |] and X(L) C It ^) '■ ~ \ — * - || contained in a 
line. 

Even though the dimension here is 2, and Lemma [3.31 applies only to dimension 1, we will still be 
able to use it for the pair (/ZB,r(L)), since the attractor X{L) is contained in a line and therefore 
Xb has only finitely many zeros in X(L) (see [DJ 07b| ). 
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Note that 



XB 



r (l + e 



2wix _|_ ^2iriy\ 



Then 



XB 



1 iff x, y £ Z, and 



XL 



XL 



(1 _|_ g27ri(a;+2j/) _|_ e -2ni(x+2y)^ 



1 iff x + 2y £ 



This implies that there are no extreme -B-cycles in X{L) and so T(L) is an ONB for fi B - Also, we 



have llial ( J is a fixed point for the map 7"(o,i) T au< l 



is a non-trivial L-extreme cycle. Hence T(B) is not an ONB for [i^. 



XL 




1/2 



1, and therefore 




1/2 



The next proposition will show that if the set L' is an integer multiple of the set L then it can 
only have more .B-extreme cycles. This implies that T(L') will have fewer chances of being an ONB. 

Proposition 3.8. Suppose (R,B,L) and (R,B,L') are as specified in Definition \2. 51 and V = qL 

for some non-zero integer q. If there are some non-trivial B-extreme cycles in X{L) then there are 
non-trivial B-extreme cycles in X(L'). Thus if the dimension is 1, and ifT{L) is not an ONB then 
r(L') is not an ONB either. 

Proof. Let C be a 2?-extreme cycle in X{L). Then, by Lemma 13.41 b ■ x S Z for all b £ B and 
x £ C. Then b ■ qx G Z for b £ B, x £ C. Therefore \xb(q^)\ = 1 for all x £ C. Also if xq,x± £ C 
and (R T )~ 1 (xo + Iq) = x\, for some Iq £ L, then (R T )~ 1 (qxo + qlo) = x\. This shows that qC is a 
i3-extreme cycle in X(L'). 

□ 



4. Transfer operators 

Recall (Corollary 12. 8p that each of the affine fractal measures under consideration has a well 
defined spectral function. Here we show (Lemma I4.ip that the pair of spectral functions corre- 
sponding to our paired fractal measures are fixed under an associated pair of transfer operators. 
Since transfer operators have a rich spectral theory, this throws light on the harmonic analysis of 
affine fractal measures. 

Let R,B,L in M rf be as specified in Definition 12.51 an d assume G B, £ L. Let (tj,)&£b and 
( T l)l&L be the dual IFSs from (|2.4|) and (|2.15p respectively. 
We will consider the following transfer operators 

(4-1) (T B f)(t) = Y,\XB(r l (t))\ 2 f(r l (t)), 

leL 

and 

(4-2) (T L f)(t) = J2\xL(n(t))\ 2 f(n(t)) 

beB 

Set 

(4.3) E + (T B ) :={f:R^ [0, 1] :T B f = f,f£ C 1 and /(0) = 1} 

and similarly for E+(Tl). 
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Lemma 4.1. Let Tg and Tl be the transfer operators in (|4.ip and (|4.2|) . Let 1 be the constant 
function one. Then 

(4.4) l€ E + (T B )nE + (T L ), 



(4.5) G E +( T b), 
and 

(4.6) a^ B) G £ + (T L ). 

Moreover, each space E + (Tb) and E + {Ti) is a convex order interval; specifically, the following 
implication holds: 

(4.7) Iffe E+(T B ) then a^ L) (t) < f(t) < 1 
and similarly for £+(Tl). 

Proof. Details are included in |DJ06bl iDTOTal IDJ07bl rPJ07cj We sketch the proof of (|I3|) . Using 
the definition of T(L) we get 

(4.8) T(L) = L + R T T(L). 
Substitution into (|2.22|) then yields 

4(i)W = E E \^B(t+i+R T ^\ 2 =j2\xB(cri(t))\ 2 E ^((^r'^+o+T)! 2 

«eL 7 Gr(L) /eL 7er(L) 



= E I^W)| 2 4(i)(^)) = 

as claimed. 

Verification of o~p,L(0) = 1: this is a direct computation using (|2.22p . and the Hadamard axiom 
in (|2.12p . See also the formula 



from Corollary! 

Indeed, setting t = into (|4.9p yields 

4f2)(°) = < e o > U (°) P LU(0)e ) L ^ B) = (e , Pzeo)^) = <eo , e ) i2(/is) = 1 

since eo is in the range of the projection P^. Recall the assumption € L. 
For the proof of gZD, see |D J06bl IDJ08] . 

□ 
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5. The Cantor measure with scale-similarity factor 4 

In this section we revisit a particular fractal measure jx. It is a Cantor measure supported by a 
compact fractal contained in the real line. Its harmonic analysis was studied first in [JP98] where 
the authors proved that it has a Fourier basis. In addition to the realization of /i as a Hutchinson 
measure, it is also (see [JP98]) the result of a recursive Cantor construction with gap-spacing and 
subdivision 4. 

Here is the specific algorithm: begin with the unit interval [0, 1], subdivide by 4, and leave two 
gaps. Then renormalize the restriction of Lebesgue measure to two of the four subintervals that 
are retained. Now continue recursively. The resulting sequence of measures has a unique limit. It 
is the Cantor measure with fractal dimension ^. 

Indeed, the limit measure [i has Hausdorff dimension ^. And the Hausdorff dimension coincides 
with the scaling dimension. The authors of |JP98j showed that L 2 (fj,) has an explicit orthonormal 
basis (ONB) of complex exponentials. The ONB found in [JP98J is here called T({0, 1}). 

In the section below we explore the possibility of scaling the earlier known ONBs by integral 
multiples p. The fact that such a scaling may even lead to new ONBs is rather surprising as the 
scaled sets become more sparse with increasing values of the integer p. Surprisingly the arithmetic 
properties of p explain and account for when the result is again an ONB. For example we show that 
if p is divisible by 3, then T({0,p}) is not an ONB in L 2 ([i). When p does not contain the prime 
factor 3, we show that T({0,p}) may or may not an ONB. For example we prove (Proposition 15.1( 1 
that the case p = 5 k is affirmative, i.e., that T({0, 5 fc }) is an ONB in L 2 {[i). 

Proposition 5.1. Let R = 4, B = {0, 2}. Then for every integer k > and L(5 k ) := {0, 5 k }, the 
setT(L(5 k )) is an ONB for L 2 (fi B ). 

Proof. Using Lemma 13.31 we will show that for any integer k > 0, there are no non-trivial B- 
extreme cycles in X(L(5 k )). We will do this by induction on k. For k = 0, 1, this is easy to check, 
see the conditions below. 

Fix k > 2. Let xq G X(L(5 k )) be a point in a I?-extreme cycle, xq ^ 0. We have 

XB (t)= l -(l + e 2 ™ 2 - t ). 

Then |xb(xo)| = 1 implies that xq G |Z. 

We claim that any such £>-extreme cycle point must be in Z. But if xq is not in Z, it must be 
of the form xq = a/2 with a odd. Let x± be the next point in the cycle, so x\ = (xo + ^o)/4 for 
some Iq £ {0,5 fe }. Since x\ is a B-extreme cycle point we also have x\ G ^Z. If Iq = 5 k then 
xi = (a + 2 ■ 5 k )/8. Since the numerator is odd, this cannot be in ^Z. If Iq = we must have 
\TL 3 x\ = | which is again impossible. Thus xq G Z. 

Let xo, a?i, . . . , x n -\ be the points in this 5-extreme cycle, xo and let Iq, . . . , l n -\ G {0, 5 fc } such 
that 

x Xq + Iq Xi + h X n _ 2 + ln-2 %n-l + h-l 

(5.1) — - — = xi, — - — = x 2 ,..., =x n _i, = x . 

Then 

xo = 4xi mod 5 fc , xi = 4x 2 mod 5 fc , . . . , x n _i = 4xo mod5 fc . 

This implies that 

(5.2) (4 n - l)xi = 0mod5 fc , (i G {0, . . . , n — 1}). 
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Let I be the largest power such that 5 divides 4 n — 1. In the case when I < k, it follows from 
(|5.2p that all elements of the cycle Xi must be divisible by 5. Then xi = 5yi for some yi G Z. 
Dividing (|5.ip by 5, it follows that yj form a I?-extreme cycle for L = {0, 5 fc_1 }, which contradicts 
the inductive hypothesis. 

So we can assume I = so 4 n — 1 is divisible by 5 fc . We will prove that in this case n is divisible 
by 5 k ~ 1 . For this it is easy to prove by induction that for any k > 0: 

(5.3) 2 4 - 5 * = l + 3-5 fc+1 mod5 fc+2 

For k = this is clear. And assuming f)5.3j) for A;, write 2 4 ' 5fc l + 3 • 5 k+1 + 5 fc+2 i then raise to the 
fifth power. Using the multinomial formula and keeping track of the terms not divisible by 5 fc+3 , 
it follows that (|5.3p is true for k + 1 as well. 

Now consider the multiplicative group (Z 5 k)* of elements in Z 5 fe that are relatively prime to 5 fc . 
Equation (|5.3p for fe — 1, shows that the order of 2 in this group divides 4 • 5 fc— 1 , and (|5.3p for — 2 
shows that the order of 2 in this group cannot divide 4 • 5 fc ~ 2 . Therefore the order of 2 is divisible 
by 5 fc_1 . 

Then, if 2 2n = lmod5 fc this implies that the order of 2 must divide In so 5 fe_1 divides n. 

But if 5 fc_1 divides n, then n > 5 . On the other hand n is the length of the cycle. This means 
that we have at least 5 fc_1 points in Z n X(L). 

However, X(L) is contained in the interval [0, 5*/4*] = [°>5 fe /3]. There are at most 5 k /3 
non-zero integers in this interval. But, as above, for all points in the cycle we must have Xi + l\ 
divisible by 4, so Xi = 0mod4 or Xi = — 5 fc mod4. These are only 2 equivalence classes mod 4, so 
the number of integers in X{L) that can be on an extreme cycle is at most (5 fc /3)/2 < 5 fc—1 . But 
since the length of the cycle , n, is at least 5 fc_1 , we reach a contradiction. 

Thus, there are no non-trivial i?-extreme cycles, and therefore T({0, 5 fc }) is an ONB. 

□ 

Remark 5.2. Proposition 15. II shows that an ONB for a fractal measure can have the corresponding 
fractional upper Beurling density arbitrarily small. 

The fractional upper Beurling density is defined as follow (see fCKS08| ): for a discrete subset A 
of M d , and for a > 0, the a-upper Beurling density of A is defined by 

X£(A) = hmsup sup — 

and the upper Beurling dimension of A is defined by 

dim + (A) = sup{a > : V a (A) > 0}. 

It was proved in [DHSW09] that, under some mild assumptions, for any set A such that the 
exponentials E(A) := {e\ : A G A} form a Bessel sequence in L 2 (u) (in particular ONBs), the 
upper Beurling dimension is equal to the Hausdorff dimension which in this case is In 2/ In 4 = 1/2, 
and the 1/2-upper Beurling density T>y 2 (A) is finite. 

Thus we have T>y 2 (Ti) < oo, where Ti := r({0, 1}). We check also that T>y 2 (Ti) > 0. 

Take x = and h = X^fc=o ^ = — -0/^ ^ n ^ e definition of the Beurling density. Then 
#(Ti n [—h, h]) = 2 n (since we have two digits in {0, 1} and n positions to write the elements in 
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IM. Then 



Thus 



P+ (r x ) > limsup /2 = > 0. 

" — X \ 3 / 



Also, it is easy to see that 



and therefore 



< V+^iTx) < oo. 

vt^) = ^^/ 2 (ri), 

lim V+ (5*Ti) = 
Let d = 1, R = 4 and B = {0, 2} so the measure fiB satisfies 

(5.4) M*) = e^ncos(^ 

k=i ^ 

Set 

(5.5) r:=r i: =|^a fe 4 fc :a fc G{0,l},nE 



.k=0 



i.e., T = T(Li) with L x = {0, 1}. 



Proposition 5.3. Let p G Z +; p > 1, and L p := {0,p}- Then the following conditions are 
equivalent: 

(i) r(Lp) = pT\ is orthogonal (not necessarily ONB) in L?(hb); 

(ii) T(B) = 2Ti is orthogonal in L 2 {^l p )- 

(iii) p is odd. 

Proof. See Proposition 12.61 □ 

Corollary 5.4. Let B and L p = {0,p} be as in Proposition 15.31 i.e., assume p is odd. Then the 
following conditions are equivalent: 

(i) The only B -extreme cycles in X{L p ) are the trivial singleton {0}. 

(ii) The only L p -extreme cycles in X{B) are the trivial singleton {0}. 

(iii) The orthogonal sets in Proposition 15.31 are ONBs. 

Proof. Combine Lemma 13.31 and Theorem 12.101 

□ 



Here is a list of non-trivial l?-extreme cycles for all the values of p < 100. This consists of all 
odd multiples of 3, and the only such p not divisible by 3 is 85. For the odd values of p that do 
not appear in the table, there are no such cycles so T(L p ) is an ONB. Therefore, the list of positive 
integers p less than 100 such that T({0,p}) is an ONB in L 2 (fi B ) is: 1, 5, 7, 11, 13, 17, 19, 23, 25, 
29, 31, 35, 37, 41, 43, 47, 49, 53, 55, 59, 61, 65, 67, 71, 73, 77, 79, 83, 89, 91, 95, 97. 
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97 
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oo 
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45 


•f 19 31 il 51 




J 1 Q 1«zl 1 \ / 1 7\ 
110,10,^,1 j-,\i ( j- 


57 






11 fi 4 1 1 11 7 2D 51- -f 21 1 


69 


{23} 


75 


{20,5},{25} 


81 


{27} 


85 


{23,27,28,7} 


87 


{29} 


93 


{31} 


99 


{33} 



Remark 5.5. It follows from the analysis of the extreme cycles introduced in Proposition 15.31 and 
15,41 that when p is divisible by 3, then there are non-trivial £>-extreme cycles in X({0,p}). And as 
a result, the orthogonal sets from Proposition 15.31 cannot be total in the respective L 2 (/i)-Hilbert 
spaces. However this implication only goes one way as is illustrated in the table for the case of 
p = 85 : Even when p does not have 3 as a prime factor there may indeed be non-trivial S-extreme 
cycles in X({0,p}). For p = 85, there is a S-extreme cycle in X({0, 85}) of length four. 

6. Finite cycles 

In this section we explore variations in the lists of i3-extreme cycles for our family of ID-examples, 
algebraic and geometric properties; and we identify two classes of such cycles. Our conclusions 
have direct relevance to orthogonal harmonic analysis of fractal measures as we developed it in 
the earlier sections of this paper. But these finite cycles are of independent interest as they have 
wider significance for factorizations of Mersenne numbers (among other topics); a subject with a 
multitude of applications: combinatorics, number theory, and encryption; see for example [FLS09, 
IMum941 iMPOl IUdl781 IVas06] . 

In addition, we mention that finite cycles that may be distributed on particular lattices have 
uses in other problems in analysis, for example in the study of representations of the Cuntz C*- 
algebras. In [BJ99], the authors introduced a family of permutative representations of each of the 
Cuntz algebras O n . They showed that finite affine cycles distributed on certain associated lattices 
account for the orthogonal decompositions of these representations. A representation of O n acting 
in a Hilbert space % is said to be permutative if it permutes the vectors in some orthonormal 
basis for %. The papers |DJ07bl DJ07c] offer yet other applications of finite cycles with lattice 
coordinates. For related questions about representations of O n , see also [SZ08J. 
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Note that in the study of fractals and dynamics, there are multitudes of families of finite cycles; 
but our present restriction that the particular cycles be .B-extreme cuts down the number of cases. 
And for the discussion below, the 5-extremality of a particular cycle C turns out to be equivalent 
to requiring C to be contained in a certain lattice. 

We consider a Hadamard system in Definition 12.51 of the following form in one dimension, i.e., 
d= 1. 

Setting. Let R = 2n, B = {0,2}, p 6 Z + odd, and set L = L„ iP = {0, np/2}. 
Then the conditions in Definition 12.51 are satisfied and the Hadamard matrix is 

-i l 1 
M 1 -h 

In this section we compute some special r^-cycles which are .B-extreme. This is of interest since 
the ONB condition in Lemma 13.31 holds iff the only .B-extreme r^-cycle is the singleton {0}. 
Recall from Lemma 13.41 that a finite r^-cycle C is .B-extreme iff C C ^Z, see also Definition 13.21 
Set 

(6.1) ToW = _ andTl(t) = _^ = _ + £ 

and let u = (uj\uj2 ■ ■ ■ ^>i), £ {0, 1} be a finite word. Set 

(6-2) T u = T U1 O T U2 O ■ ■ ■ O T Ul . 

Lemma 6.1. If p is divisible by 2n — 1 then there are non-trivial B-extreme r^ n p -cycles of length 
one. 

Proof. Let p = m(2n — 1), m E Z. Then the solution t to T\{t) = t is 

np mn 1 „ 

(6.3 t = — — — -- = G -Z. 

v ; 2(2n - 1) 2 2 

Hence {^} is a l?-extreme tl 71 p -cycle. □ 

Corollary 6.2. (i) For R = 4, B = {0,2}, and L = {0,3}, the singleton C = {1} is a B-extreme 
TL-cycle; and so in particular T({0,3}) is not an ONB in L 2 (// 1 / 4 ). 

(ii) For R = 6, B = {0, 2} and L = |0, the singleton C = {|} is a B-extreme TL-cycle. 

(Hi) For R = 8, B = {0, 2}, and L = {0, 14}, the singleton C = {2} is a B-extreme TL-cycle. 

Definition 6.3. Let n £ Z + be given, and set R = 2n, B = {0,2} and let fi = ^i/2n = Mb be 
the fractal measure from section O We say that p G Z + odd is admissible iff there are non-trivial 
.B-extreme ?x-cycles. Here L = L n ^ p = {0, np/2}. 

Note that if n = 2 then p = 3 is admissible. A consequence of Lemma 16.11 is that for n G Z+, 
p = 2n — 1 is admissible, and, with Proposition I3.8( we obtain that any p divisible by 2n — 1 is 
admissible. 

Question. Are there admissible values of p not divisible by 2n — 1? 
The next theorem offers an affirmative answer. 

Theorem 6.4. Let n G Z + be given. Set 

2n-l 

(6.4) p=£(2n) < . 

i=0 
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Then p is admissible and not divisible by 2n — 1. There are associated B -extreme cycles of length 
2n. 



Proof. Let p be given by (|6.4p . Then p = lmod2n — 1, so it is not divisible by 2n — 1. We shall 
need 

2n-2 

(6.5) p* := ^ (2n)\ 

i=0 

Note that p* = 0mod2n — 1. So p* = m(2n — 1) for some m£Z. 
Consider the solution t of the fixed point equation 

It is easy to see that 

_ IT X^=o~ 2 ( 2n )' _ npp* np* _ nm 1 

(2n) 2n - 1 2(2n - l)p 2(2n - 1) 2 2 

We check by induction that (t) G |Z for fc G {0, 1, . . . , 2n - 1}. We use that p = p* + (2n) 2n ~ 1 . 
We have 

r(t) _t + ! ¥_ T + f = n ^ =l +V* + {2nf^ = n ( (2n)p* 2n _ 
U ' 2n 2n 2 2n 2 \2n{2n - I) K ' 



So 

Then 



n{t) = — + -{2nY n ^ = t + -{2n 



n >~ ^2n-2 



\t) = Tl (t + ^(2n) 2 - 2 ) = Tl (t) + ^^(2n) 2 - 2 = t + \(2nf^ + ^(2n)^ 



where we used the previous step in the last equality. 
By induction we get 

Tl fc (t) = t + \{2nf n - 2 + \{2nf^ + ■■■ + ^nf^ 1 

for all k < 2n — 1. And this shows that r^(i) is in iZ. Since ror 2n_1 (t) = t, this implies that the 
entire cycle is in |Z, so it is a 5-extreme Ti np -cycle. 

□ 

Remark 6.5. An application of Theorem 16.41 to R = 4 and p = 1 + 4 + 4 2 + 4 3 = 85 accounts for 
the cycle C = {23, 27, 28, 7} from the table before Remark 15.51 It is the smallest admissible p not 
divisible by 3. 

For R = 6 and p = 1 + 6 + 6 2 + 6 3 + 6 4 + 6 5 = 9331 on has L = {0, 27993/2} the 5-extreme cycle 
of length 6 is 



f 4821 5469 5577 5595 9331 

C = { , , , ,2799, } 

\ 2 ' 2 ' 2 ' 2 ' ' 2 J 



It is also the smallest admissible p not divisible by 5. 

For R = 8 and p = ^J=q 8* = 2396745 one has L = {0, 4793490} and the ^-extreme cycle of 
length 8 is 

C = {609886, 675422, 683614, 684638, 684766, 684782, 684784, 85598} . 
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We used the following Mathematica program to check for extreme cycles, 
r = 8; 

NextCycle [x_ , p_] : = 

If [IntegerQ [2*x/r] , x/r, If [IntegerQ [2* (x + p)/r] , (x + p)/r, -1]]; 
pmax = r*(r~r - l)/(4*(r - 1)) 
For[p = pmax - r/2, p < pmax, p = p + r/2; 
If[! IntegerQ [2*p/(r - 1)], 
For[ix = 1/2, ix <= IntegerPart [2*p/ (r - l)]/2, ix = ix + 1/2, 
x = NextCycle [ix, p] ; flag = 0; 
While[x != -1 && x != ix, x = NextCycle[x, p]]; 

If[x != -1, cycle = {}; If [flag == 0, Print ["p= ", p] ; flag = 1]; 
Print ["Cycle pt "]; i = 1; 

While[i == 1, x = NextCycle[x, p] ; Print[" ", x] ; 
If[x == ix, i = 0]]; Print ["End Cycle "] 

] 

] 

] 

] 
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